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PRECISION ADAPTABLE CHEMICAL SPECTROMETER 
Reference to Related Applications 

This application claims priority to U.S. Provisional Application No. 60/391,770, 
entitled "Integrated Remote Chemical Spectrometer" filed June 27, 200, which is hereby 
5 entirely and completely incorporated by reference. 
Background 

1. Field of the Invention 

This invention is directed to spectrometers that fosters and detects absorptions in 
submillimeter wave region, and to methods for processing spectroscopic data with unique 
10 algorithm. 

2. Description of the Background 

Spectroscopy using submillimeter wave, or microwave, has the potential to be a 
potent diagnostic tool. Molecules having a permanent or induced dipole moment have 
discrete rotational energy levels that can absorb electromagnetic waves in MHz to THz 

15 regions. Because of the abundance of rotational energy levels that are unique to the 
structure of the molecule detected, submillimeter spectroscopy makes it possible to 
identify chemical species with excellent specificity. One of the earliest studies in this 
filed was reported almost 50 years ago; "The use of microwave spectroscopy for 
chemical analysis," C. H. Townes and A. L. Schawlow, chapter 18, "Microwave 

20 Spectroscopy," McGraw-Hill, 1995. However, this analytical method did not result in a 
wide spread use in research communities and industries, partially because the complexity 
of the instrumentation and associated high cost and maintenance issues. 

Recently, with the advancement of THz technologies and computer based system 
control, submillimeter spectrometers with more practical use have been developed. Such 

25 spectrometers are reported, for example, in "A Fast Scan Submillimeter Spectroscopic 
Technique," D.T. Petkie, T.M. Goyette, R.P.A. Bettens, S.P. Belov, S. Albert, P. 
Helminger, and F.C. De Lucia, Rev. Scient. Instrum 68, 1675-1683 (1997), and 
"FASSST: A new Gas-Phase Analytical Tool," S. Albert, D.T. Petkie, R.P.A. Bettens, 
S.P. Belov, and F.C. De Lucia, Anal. Chem. 70, 719A-727A (1998). 

30 This type of spectroscopic technique was initially made possible by a combination 

of fast scanning voltage tunable Backward Wave Oscillators (BWOs), optical calibration 
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methods, and modern fast digitization and computation techniques. Fig. 1 shows a 
diagram of such a spectrometer. Briefly, the BWOs are voltage tunable (-1500 - 4000 
V) tubes covering roughly 0.1 - 1 .0 THz, in bands. A typical tube, an OB-30, covers 
-250 - 350 GHz as the voltage is swept over this range. In this system the output of the 
5 BWO is split, with about 10% going to a Fabry-Perot cavity to provide optical 

calibration, while the rest is used to interrogate the sample. The portion of the power 
then passes through the sample cell and is detected by an InSb detector. In operation, the 
tuning voltage is ramped to provide an analog sweep in frequency, and the outputs of the 
signal detector and FP detector are recorded in parallel. 

10 A typical sweep time is one second and the frequency width 10-100 GHz. This 

frequency interval contains ~10 5 independent resolution elements. With integration times 
of 1 microsecond, S/N ratios of 10 4 are obtained. 

The system concept is based on the short term spectral purity of the BWO (Q > 
10 7 ). Because of this spectral purity, the slow and complex phase lock that is ordinarily 

1 5 used to stabilize and control THz spectroscopic systems is not fundamentally necessary. 
In this spectrometer, the high spectral purity source is swept so rapidly that slow 
instabilities (associated with thermal drift, etc) are frozen on the time scale of a 
measurement cycle. Fast digitization records the output of the FP cavity and 
spectrometer in parallel and makes possible accurate frequency calibration even with 

20 thermal drifts, power supply ripple, and nonlinear frequency sweeps. In some sense, the 
speed of the digitization plays the same role as the bandwidth of a more traditional phase 
lock loop. 

The scanning speed combined with the high spectral brightness (W/Hz) of 
electronic sources and the very strong molecular interactions in the THz range, rapidly 
25 produces analytical fingerprints with remarkable information content. This is illustrated 
in Fig. 2, which is presented as a series of blow-ups in both frequency and amplitude. 
Fig. 2 shows about a 25 GHz region of the spectrum of HN03, which was recorded in -1 
second. The middle and lower figures show expansions of this spectrum. 

Stated another way, if the upper figure were plotted with 1 mm of noise in the 
30 vertical and 1 cm of width per resolution element in the horizontal, it would require a 
piece of paper 10 m high and 1 km long for each second of data acquisition. 
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The high resolution of the millimeter spectroscopy comes about because the 
Doppler broadening (which is the fundamental limit in most systems) is proportional to 
frequency. Thus, lines in THz spectra are about 100 times narrower than those in 
infrared spectra. For example, for molecules of about the size of CIONO2, their spectra 
5 are unresolvable in the infrared. Rather than the detailed rotational spectra shown in Fig. 
3, only a broad band is observable in a infrared spectrum. Because of the high resolution 
in the THz, the optimum pressure (set by the condition that the Doppler and pressure 
broadening be approximately equal) is much lower, typically 10-50 mTorr. 

One of the disadvantages of this spectrometer is, however, its size. The Fabry- 

10 Perot cavity alone occupies a space as large as a small room. In addition, BWO and 
associated magnets for guiding electrons as well as the InSb electron bolometer that 
needs to operate at 1 .5 K require plenty of space for its instrumentation. Because of its 
excellent specificity and speed, the submillimeter spectrometer is theoretically suited for 
applications outside laboratories such as remote chemical analysis. However, at least the 

1 5 large size has prevented the submillimeter spectrometer from being seriously considered 
for such applications. 
Summary of the Invention 

The invention overcomes the problems and disadvantages associated with current 
designs of the submillimeter spectrometers and provides a new compact submillimeter 

20 spectrometers and methods for identifying chemical species using the spectrometer. 

One embodiment of the invention is directed to a spectrometer for measuring 
submillimeter absorption. The spectrometer includes a solid state exciter generating a 
submillimeter wave and sweeping a predetermined frequency band, a frequency marker 
generating unit electrically generating frequency markers, a sample cell to contain a gas, 

25 and a solid state detector detecting a submillimeter absorption of the gas. 

Another embodiment of the invention is directed to a spectroscopic method that 
includes energizing a solid sate oscillator to generate a submillimeter wave and to sweep 
a predetermined frequency band, introducing the submillimeter wave during the sweep 
into a sample cell that contains a gas, electrically generating frequency markers 

30 sequentially during the sweep, reading during the sweep outputs of a solid state detector 
disposed in the sample cell and detecting an absorption of the gas, recording the read 
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outputs of the solid state detector as a function of time and recording the frequency 
markers as a function of time, and converting the recorded outputs of the solid state 
detector into a function of frequency using the recorded frequency markers. 

Another embodiment of the invention is directed to a method of measuring 
5 submillimeter absorption. The method includes introducing a gas containing a chemical 
species of interest into a trap, cooling the trap so that the chemical species is in a 
condensed state, reduce the pressure of the trap to evacuate the remainder of the gas that 
is not condensed, heating the trap to a first temperature so that the chemical species 
evaporates, leading the evaporated chemical species into a sample cell, and measuring 

1 0 absorption of the chemical species. 

Another embodiment of the invention is directed to a method of surveying an area 
for a chemical species. The method includes taking air into a spectrometer at a first 
location, receiving global positioning system (GPS) coordinates at the first location, 
measuring a submillimeter absorption spectrum of the air taken in at the first location, 

1 5 recording the absorption spectrum of the first location with the GPS coordinates of the 
first location, taking air into a spectrometer at a second location, receiving GPS 
coordinates at the second location, measuring a submillimeter absorption spectrum of the 
air taken in at the second location, and recording the absorption spectrum of the second 
location with the GPS coordinates of the second location. 

20 Yet, another embodiment of the invention is directed to a method of determining a 

presence of a chemical species in a mixture of gasses. The method includes introducing 
the mixture of gasses into a sample cell, irradiating the mixture of gasses in the sample 
cell with a submillimeter wave that sweeps a predetermined frequency band, generating a 
submillimeter spectrum of the mixture of gasses, providing a standard submillimeter 

25 spectrum of the chemical species that is obtained from the chemical species of a pure 
form, selecting a first peak of the standard submillimeter spectrum of the chemical 
species, and determining whether the selected first peak is present in the generated 
submillimeter spectrum of the mixture of gasses. 

Other embodiments and advantages of the invention are set forth in the 

30 description, which follows, and in part, may be obvious from this description, or may be 
learned from the practice of the invention. 
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Description of the Figures 

Figure 1 is a system diagram of a conventional submillimeter spectrometer. 
Figure 2 is an example of spectra taken with the spectrometer of Fig. 1. 
Figure 3 is another example of the spectra taken with the spectrometer of Fig. 1 . 
5 Figure 4 compares spectra of a gas mixture taken with and without a specific chemical 
species. 

Figure 5 is a system diagram of a submillimeter spectrometer of an embodiment of this 
invention. 

Figure 6 shows the frequency marker generation unit of the system shown in Fig. 5. 
10 Figure 7 is a schematic diagram of a gas processing system of the spectrometer of Fig. 5. 

Figure 8 is a schematic drawing to show the effect of the removal of ambient gas 

components on the spectrum sensitivity. 

Figure 9 schematically show a fractionation loop. 

Figure 10 is a diagram of the signal processing scheme. 
1 5 Figure 1 1 shows the probability of detection as a function of probability of false alarm. 

Description of the Invention 

As embodied and broadly described herein, the invention is directed to a method 

for processing spectroscopic data with unique algorithms and a compact submillimeter 

spectrometer based on these algorithms. 
20 A submillimeter spectrometer, such as shown in Fig. 1, has been used in 

laboratory settings to investigate absorptions in submillimeter wave region. That is, a 

fluid sample, which may be a gas, a liquid, or an aerosol, is prepared prior to the 

spectroscopic measurement so that the fluid contains only a single chemical species or 

majority of the fluid is the single chemical species. However, the submillimeter 
25 spectrometer has not previously been utilized to identify a single chemical spices in a 

fluid mixture that contains the chemical species as well as other major and minor 

ingredients. 

Fig. 4 is a result of spectroscopic measurements directed to identifying a single 
chemical species in a gas mixture. A segmented sample cell with twenty-one (21) 
30 separate vacuum segments along the optical path, was constructed to keep separate the 
various possibly reactive chemical species within a sample. This sample cell was 
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incorporated into the submillimeter spectrometer shown in Fig 1 . Cells may further 
contain one or more agents to inhibit interaction between chemical species with a sample. 
Agents may be incorporated into the cell structure, layered onto the inner surface of the 
cell, or added to the sample at any stage. Agents may be general inhibitors of all 
5 chemical interaction, or specifically targeted to inhibit only specific chemical reactions. 
Alternatively, the sample within the cell may be exposed to conditions to minimize or 
substantially inhibit general or specific interaction between chemical species of interest. 
Such conditions include, but are not limited to, extremes in temperature, pressure or 
exposure to radiation of defined wavelengths. 

1 0 In this experiment, 30 chemically library spectra were provided and compared to 

composite spectra of 19 chemical spices. In addition, the first segment was used to 
contain SO2 for registering the Fabry-Perot fringes, and the last segment contained the 
atmospheric "soup" for the clutter aspect of the experiment. The target chemical spices 
of the experiments shown in Fig. 4 were CH2F2 and vinyl chloride. Despite the presence 

15 of other 18 chemical spices in the gas, the presence or absence of the species of interest 
from the overall mixture was confirmed. There was an immense degree of spectral 
information content in the top overall spectrum for each gas. The lower zoomed views, 
which show only about one percent (1%) of the data on a still compressed scale, provided 
easy spectral identification substantially without influence of other species present. 

20 Many such intervals can be examined with similar conclusions. The information is very 
abundant so that it is essentially multiply redundant. The abundance of usable spectral 
information provides redundancy in detection/identification opportunities, thus providing 
high sensor performance in false alarms (i.e., misidentifications) and detection 
probability (i.e., given the true presence of a candidate species). This spectrometer 

25 system may be used to create a spectrum library that is used with a compact spectrometer, 
such as described herein. 

Fig. 5 shows the overall system design of a submillimeter spectrometer of an 
embodiment of the invention. The system is divided into two subsystems, i.e., THz 
module subsystem and digital processing and control subsystem. The THz subsystem 

30 includes a gas process unit having a gas trap, a sample cell for spectroscopic 

measurement and gas outlet unit as well as a solid state exciter coupled with an electronic 
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reference system. A gas to be examined, or most likely a mixture of gasses that may 
include a gas of interest, is led into the gas trap. The gas is trapped in the gas trap under 
cooling to separate it from other chemical elements of the mixture, and then heated to 
evaporate into the sample cell. The all solid-state exciter generates a submillimeter 
5 incident wave of 300 GHz, and, at the same time, reference markers for calibration are 
generated by the associated electronic reference system. The frequency markers are 
generally of the same frequency throughout a measurement. Exact value of the frequency 
marker may depend on the overall system design. For example, it may be 100 KHz, 1 
MHz, 10 MHz, 25 MHz, 35 MHz or the like. The signals from a diode detector located 

10 in the sample cell and the corresponding reference markers are stored in the digital 
processing and control subsystem. This subsystem also controls the exciter, the 
electronic reference system and the gas process unit. 

The principle of detecting absorption peaks corresponding to the rotational energy 
level(s) is the same as that of the spectrometer with the B WO and Fabry-Perot cavity 

1 5 shown in Fig. 1 . However, in this embodiment, the radiation source, B WO, is replaced 
by the compact solid state exciter, and the optical reference system, Fabry-Perot cavity, is 
replaced by the electronic reference system. In addition, the InSb electron bolometer is 
also replaced by the Schottky diode. Accordingly, the overall size of the spectrometer is 
significantly reduced. Preferably, the spectrometer of this invention has a size small 

20 enough to be deployed in a wide variety of field applications and may be portable such 
that it can be carried by a single individual, and easily packaged. Suitable sizes are about 
five (5) square meters, about two (2) square meters, about one (1) square meter, about 
one half (1/2) a square meter, about one hundred (100) square centimeters, or even less. 
Suitable weights are less than two hundred (200) kilograms, less than one hundred (100) 

25 kilograms, less than fifty (50) kilograms, less than twenty five (25) kilograms, less than 
ten (10) kilograms, and less than five (5) kilograms. 

Components of the solid-state exciter may be chosen from those currently 
available commercially. The generation of the radiation wave starts with a 10-20 GHz 
VCOora5-10GHzVCO. In addition VCXOs may be used. Several VCOs were tested 

30 and found to have the required spectral purity (<10~ 6 ) and tuning range. These VCOs 
have considerably more linear tuning than the BWOs, which simplifies the frequency 
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calibration. The VCO of 1 mW drives an active multiplier chain to produce -40 mW 
near 100 GHz, followed either by one stage of a passive diode frequency doubler or 
tripler, or two frequency doublers in tandem. This system will occupy a few 10s of cm 3 
and require ~10V. 

5 With 1 mW of available source power, the InSb bolometer detectors are driven far 

into saturation. Schottky diode detectors are available with responsivities of > 1000 V/W. 
With the drive power available, performances are similar to those of helium detectors. A 
less sophisticated diode detector at ~1 50 GHz was tested and still found to be superior to 
the cooled bolometer detectors because it saturates at a higher power level. The Schottky 
10 diode shown in Fig. 5 has 1000 mV/mW. Because of the use of solid-state devices in 
place of Fabry-Perot cavity, BWO and the bolometer, the whole spectrometer system 
occupies a volume on the order of <1 ft 3 with power consumption of <100 W for the 
entire system. However, a spectrometer with a volume less than 5 ft 3 or 10 ft 3 would also 
be useful. 

1 5 The operation of the THz module is further described herein. The controller of 

the digital processing and control subsystem sends an approximately linear ramp to the 
10 GHz VCO of the THz module, which responds with an approximately linear analog 
sweep of frequency as a function of time, typically taking ~ 1 second to produce spectra 
similar to that shown in top of Fig. 2 or -0.1 and -0.01 seconds for the smaller subscans 

20 shown in Fig. 2. In the example shown in Fig. 5, three multipliers work together to raise 
the frequency of the incident wave to about 300 GHz for a sweep of a 30 GHz band. The 
generated submillimeter wave is guided to one end of the sample cell, partially absorbed 
by the gas species while traveling through the cell that has an absorption effective path 
length of about one (1) m, and detected by the Schottky diode at the other end of the cell. 

25 The frequency markers are generated based on the output of the same 10 GHz oscillator. 
The integration time for the signal acquisition is about 1 microsecond in this example. In 
this example, references markers are sequentially generated during the sweep of the 30 
GHz band. This corresponds to about 30,000 resolution elements at 300 GHz. The 
controller digitizes both the output of the spectral channel and frequency marker channel 

30 in parallel, and software uses the information in the marker channel to assign a frequency 
to each time/frequency bin in the spectrum. Thus, simultaneous collection of reference 
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markers within the 30 GHz band pass of the spectrometer allows the time scale to be 
related back to frequency. This spectrum, calibrated in both frequency and intensity, can 
then be compared manually with a library, or be used in conjunction with sophisticated 
algorithms for probability of detection determinations. 
5 Alternatively, there may be a particular gas or gasses of interest and it may be 

advantageous for the system to spend observation time focused on that gas. In this case, 
the system can focus its frequency sweep on a subscan or subscans of the entire spectral 
interval to either achieve greater speed or sensitivity. Because even these subscans 
contain many resolution elements, in most cases false alarm rates of near zero are still 
10 obtained. This is useful especially when the spectrometer system is not able to isolate a 
chemical species from the mixture taken into the spectrometer. In such an application, a 
frequency subband is positioned to avoid absorption peaks due to other ingredients of the 
gas mixture. The width of the subband is, for example, 10 GHz, 5 GHz, 1 GHz, 100 KHz 
or the like. 

1 5 The acquired data produces lines at essentially their Doppler line widths (about 1 

MHz). This spectrum is digitized and stored in memory at a rate of —106 data points/sec. 
Because the optimum pressure for spectroscopy in the THz is low (-10 rnTorr or about 
10* 5 atmospheres) and because the eventual sample cell volume can be made as small as 
-1 cm 3 , the total sample required is very small. If this small sample is not diluted in a 

20 background atmosphere, the signals are very strong, even with the 10" 6 second integrating 
time. 

Fig. 6 shows further details of the frequency marker generation scheme shown in 
Fig. 5. In this example, the scanning VCO is mixed with a crystal based reference at -10 
GHz. The mixing products (± ~1 GHz) is amplified and mixed with a 10 MHz comb. 

25 The mixing product (±100 kHz) is then amplified and digitized to produce absolute 
frequency reference points, by interpolation, for the entire spectrum. In addition to the 
very large size reduction, this is attractive because the solid state system starts with a 
VCO at a relatively low frequency (-10 GHz). This allows the use of standard 
microwave components to obtain a frequency reference. However, because of the fast 

30 sweep rate, phase information will need to be recovered and processed to provide 
accurate frequency reference. 
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Operation of the spectrometer is described with respect to the digital processing 
and control subsystem. Running a scheduler, the master controller (MC) autonomously 
regulates the gas processing system to collect a sample starting at a scheduled time or 
GPS location. The MC opens the collection and outlet valves and enables the cooler 
5 within the thermal trap while storing the GPS coordinates. Low vapor pressure samples 
condense in the thermal trap while the volatiles such as oxygen and nitrogen are pumped 
through and out of the system. The MC then closes the sample inlet and proceeds to 
pump down the remaining gas in the sample cell. Upon closing the outlet valve the 
thermal trap can be heated to boil off all or some of the sample material. Through closed 

1 0 loop monitoring of temperature and pressure data, heat can be applied to maintain a 

constant vapor pressure of the sample being measured. In this way the MC controls data 
collection through precision distillation of the collected sample. This allows the 
collection of several data sets from the same sample at varying vapor pressures. 

Once the MC has prepared the sample cell, the MC issues a command to the 

1 5 digital signal processor (DSP) to start data collection. Data collection is initiated when 
the DSP writes a new output level to the DA converter. This output level parameter is 
passed to the DSP from the MC in the "start collection" command along with the pulse 
duration. The change in voltage at the output of the AD causes an integrator in the solid- 
state exciter to start ramping the VCO, which generates the instrument frequency sweep. 

20 At the same time the frequency sweep is induced, the DSP starts simultaneous data 

collection of the detector voltage and the reference markers. Both channels can be read 
one million times over the duration of the sweep as defined in the "start collection" 
command. At the end of the frequency sweep, data collection stops and the output of the 
DA can be set low. The DSP now processes the data with several options possible. 

25 First, the DSP processes the reference data to locate the reference marker peaks 

and tags the corresponding detector output. This reduces the data set by fifty percent 
(50%) while maintaining frequency information. Since the 12 bit data sample takes up a 
byte and a nibble, the remaining nibble can be used to encode the frequency information 
directly into the detector sample. The data is then transferred to the MC where it is 

30 wrapped with header information containing GPS location and time. 
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Second, the DSP processes the reference markers as described in first scenario 
and tags the corresponding detector output. Additionally, the DSP finds the top 30,000 
peaks (or fewer if user defined and there are likely to be only a few thousand peaks in the 
30 GHz sample) and reports back offset and amplitude sets for each value. The offset is 
5 the number of whole marker ticks from the start of the sample along with the fractional 
tick information in the way of the number of samples past the tick. This reduces the data 
set to Jess than five percent (5%) of the first scenario. Note that the maximum number of 
resolution elements in a 30 GHz region around 300 GHz is about 30,000. This is a 
completely filled spectrum. The data rate is higher because it takes about ten (10) 
10 samples/line width to get some idea of line shape and to help separating signal from 
noise. Thus, the approach reduces the amount of data by at least a factor of ten (10) or 
more. 

Third, the DSP processes the samples and corrects for frequency sweep non- 
linearity as described. The processed sample is then compared to a reference library of 
1 5 compounds. Correlation of the current sample to a stored reference defines an alarm 
condition and the compound information is transmitted to an examination station. 

The digital processing and control subsystem is designed to accommodate 
operation in any one of the above scenarios, providing design flexibility. 

Fig. 7 shows an example of the gas processing unit shown in Fig. 5. VAR speed 
20 fan or a forepump is used to take the gas samples into the thermal trap, and micro 

thermoelectric coolers are used to freeze the gas samples in the trap. The optical path for 
the spectroscopy is formed as a waveguide and has an effective optical path of about 1 m 
and 1 cm 2 cross section in the sample cell. After the measurement, the gas is evacuated 
to the outside by a micro vacuum pump. In addition, the detectors and the exciters may 
25 be as small as 0.5 cm in length. 

In most applications, the gases of interest are diluted in air and are initially 
ingested. While it is still possible to detect the diluted sample by simply reducing the 
total pressure to -10 mTorr (a pressure that is obtainable with simple pumping systems 
because of the small cell volumes), it is advantageous to selectively remove the 
30 components of the ambient air which are not of interest. While a number of collection 
and concentration approaches have been developed, thermal fractionation used in the gas 
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processing unit shown in Fig. 7 is both a practical and conceptually straightforward 
approach. Air is -99% N2 and 02, and since these gases are both much more volatile 
than virtually all of the gases of interest, a simple trapping and distilling cycle can 
eliminate the air. Fig. 8 schematically shows the effect of the removal of the ambient 
5 gasses on the submillimeter spectrum. During the cool down and the associated pressure 
reduction, baseline noise level was clearly reduced. In addition, some of the peaks 
present at the beginning of the cool down disappeared at the end of the process. 

While this is a routine operation in the laboratory and even in some field 
instruments, the spectrometer of this embodiment uses more progressive fractionation to 

10 obtain a fast cycle time so as to produce as small a sampling cycle time as possible so that 
many samples can be taken quickly. This approach is especially attractive when only 
small samples are requited. 

This sampling cycle is shown in Fig. 9. The processing may best be considered as 
a series of "fraction" cycles (the lower loop shown in Fig. 9) operating within a larger 

1 5 loop denoted as the integration and reporting cycle in the same illustration. The gas 

intake system opens up and pumping begins as the thermal trap is cooled down. When an 
adequate sample is ingested and the gases of interest are frozen out, the gas inlet valve is 
closed and pumping removes the nitrogen and oxygen, which if present, induces 
undesirable pressure broadening (or alternatively, dilution if the total pressure is simply 

20 reduced by pumping) of the spectral lines as shown in Fig. 8. Next, the gas exit valve is 
closed and the frozen sample is warmed slightly (first stage "fraction") and some gas is 
created that is then probed with radiation from the solid state emitter, thus forming a 
spectral data set. With the formation of each spectrum, the markers are utilized to 
interpolatively form a uniform time/frequency baseline. This also facilitates multi scan 

25 integration. The cell is flushed and closed and the steps repeated, incrementally 
increasing the temperature and producing separate spectra for each set of molecular 
species contained in the respective "fraction." When these cycles are completed, the cell 
is purged, the data is transmitted to the analyzer and the system prepares to ingest another 
atmospheric sample. Depending upon the volatility of the gases of interest, the range of 

30 temperatures in this cycle can be adjusted in software. In addition, the "fractions" are 
selected in terms of freezing points (boiling points) of each species in the chemical 
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library associated with the spectrometer, or other criteria based on specific requirements 
of the application. For most applications, the temperatures can be obtained with small 
thermoelectric coolers. While this is conceptually the simplest gas capture front end for 
the spectrometer, its modular design makes the integration of other approaches (including 
5 those supplied by end users) straightforward. In remote mobile applications, the resulting 
composite sample can be tagged with GPS position and time based on platform 
navigation input or even an internal GPS receiver. The simplified loops in Fig. 9 do not 
show the integration of multiple scans within a particular fraction, however, this is a 
sensitivity enhancing option. Sensitivity of the submillimeter spectrometer may be 

10 further enhanced by conventional approaches as well as approaches described herein. 

First, while there are applications (e.g., the analysis of chemicals from a vessel, 
munitions, etc.) that produce a concentrated sample of the analyte, there are many 
important applications in which the chemicals are dispersed in air. In these cases the 
ultimate sensitivity of any gas analysis device can be increased if the major constituents 

1 5 of air are first removed as described herein. This approach is especially attractive in the 
THz, because the optimum sample pressure is ~10" 5 atmospheres, thus leading to 
proportionally small sample requirements. 

The major constituents of air, nitrogen and oxygen, are relatively easy to remove 
because their high volatility and small physical size are rather different from gas samples 

20 of interest. If the air in the sample cell is drawn through the trap whose temperature is 
>77K, the oxygen and nitrogen will pass and the remaining fraction will be trapped. The 
pumping valve can then be closed and the trapping surface warmed. In the simplest 
scheme the trap is warmed to 300 K and all of the trapped sample released into the cell. 
Depending upon cell volumes, pumping speeds, etc. this cycle can be done on a one (1) 

25 second to one (1) minute time scale. There are reports in the literature of this thermal 
cycling being done on time scales of a few milliseconds. 

Even in this simple example, the sensitivity is increased by a factor of 100. Only 
a modestly more complex algorithm makes it possible to look at the fractions that come 
out of the trap as a function of its temperature as described. While there are limitations 

30 on the separation of gases with similar vapor pressures, there are very large, but scenario 
dependent, gains in detection that can be achieved. 
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Second, the most straightforward enhancement of the sensitivity is to increase the 
integration time from one (1) microsecond to a longer period. There are many 
applications for which it will be unnecessary to look at all 30,000 channels to get absolute 
specificity and generality. In these, it would be advantageous to look at only a few 
5 channels, but with longer integration times. This gain may ultimately be limited by 
systematic effects that can not be reduced by averaging the spectra. Experience with 
similar systems has shown that the principle systematic effect comes from baseline 
variations. These can be traced to multipath interference between different modes as 
power is transmitted from the THz source to the detector. At least one (1) second 

1 0 integration times is expected to be useful. 

There is an even more important issue, optimizing at one (1) sec/resolution 
element, the trade-offs between speed, sensitivity, and specificity. If one were to 
integrate the entire spectrum, it would require 106 seconds. However, this is not really 
necessary and much more optimum approaches exist. For example, most of the weaker 

1 5 species are weak because their spectral intensity is divided among many lines. The 
automatic recognition software may, in effect, "integrate" these lines together. 
Accordingly, the long signal integration time for enhancing the sensitivity may not be 
necessary. 

Third, while the one (1) meter effective path length (EPL) of the sample cell is 
20 similar to those which have been used for many years, there are enhancements possible to 
overall system performance that utilize in one of several fashions multi-pass cells to 
lengthen the EPL for the absorption. The sensitivity of the system is proportional to the 
path length of the absorption cell. While this cell can be very simple (e.g., a hollow 
metal pipe 5 mm in diameter), it would be advantageous to have the absorption length to 
25 be as long as possible. This approach may result in an enhancement by a factor often 
(10), one hundred (100) or more. 

Fourth, it is useful to optimize the noise spectrum of the room temperature diode 
detectors and the preamplifier frequency response relative to the system sweep speed 
(and thus its central data recovery frequency). At a fundamental level a gain here is 
30 possible, but unpredictable system noise is an issue and accordingly the potential 

enhancement for this is estimated as ten (10). Additionally, with the availability of more 
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power in solid state THz exciters, 'autodyne' detection techniques may be used to increase 
detection sensitivity. Briefly put, power that is greater than that required for molecular 
saturation can be separately transmitted to the detector, which then plays the role of a 
mixer (with zero-frequency IF) in a heterodyne system. A similar scheme has been used 
5 for many years at lower microwave frequencies in radar speed detection devices. 

The signal processing of the spectrometer is further described with reference to 
Fig. 10. The absorption spectrum is formed through the square law response of the diode 
detector. The fluctuations in voltage from this device are very rapid by virtue of the 
sweep rate (e.g., 30 GHz/second). The system instabilities are lessened and additionally 

10 a high pass (mid frequency range) filter may be used with cut-off above the practical 
tuning margin. This yields wobbling baseline reduction and forms first derivative 
response, which is desirous. Within the detection block, residual RF is cleared prior to 
application of an LNA. The signal is then digitized through an AD converter. In parallel, 
the harmonic frequency marker generation circuit produces bit markers in the output that 

1 5 tag the exact positions of a comb line of 1 0 MHz separations. To complete this stage 
(segment A of Fig. 10) the software resolves non-linearities in the sweep, interpolating 
the equi-frequency markers to smooth the now digitized envelope. Implicitly, in this 
section of the process, built-in spectrometer calibration curves are used to adjust the 
signal amplitude in order to facilitate eventual concentration measurements, and an 

20 absolute frequency anchor is assigned to the comb of marker frequencies. This 
compensates for non-uniform amplitudes presumed here to be a function only of 
frequency. These activities are also needed at this processing level to facilitate 
integration, shown in the figure. 

As shown in segment B of Fig. 10, the system then operates upon the waveform 

25 of first derivative spectral lines by convoluting a unit Lorentzian first derivative with the 
spectrum using a limited extent of the integral with the extent determined by the level of 
smoothing desired and the processing time consumed. The resulting lines are essentially 
second derivatives, the smooth peaks of which are recorded (f,A). A clutter avoidance 
routine executes next that considers the local distortion of the spectrum at a particular line 

30 position by the presence of other spectral lines. The level of disturbance (and thus 

diminished reliability in amplitude measurements and to an extent also frequency) can be 
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determined and used to define the clutter avoidance radius that in turn operates the clutter 
avoidance logic. Lines may be denoted as "disturbed" and that factor is taken into 
account in subsequent processing. By the end of this process, the discreet spectrum is 
formed: {fk, Ak}. The final portion of the chain is shown in segment C of Fig. 10. Here, 
5 whatever operative library is in effect is summoned and run against the discreet spectrum. 
This processing is serial and includes an important clutter abatement feature. Note that as 
in the initial RF sweep, the swept or examined bands are set by operational programming. 
In this case the library contains recommended sub-bands to utilize. This will reduce 
computation loads. 

10 The ultimate assessment stage may be performed in an external PC, but is 

preferably performed internally in the spectrometer. The baseline algorithm is developed 
specifically to work in adverse conditions, i.e., many clutter lines and large dynamic 
range between various lines. The algorithm employs the unique pattern of the library 
spectrum to detect and verify that pattern within the discreet spectrum. In this process, 

15 interference by clutter lines is accounted for and generally defeated, with degradation 
only when the number of effective usable library lines is too far reduced. Eventually, 
even this degradation will normally be removed by additional subband examinations. 
The procedure is piecewise integrable so that there is absolutely no penalty in employing 
any number of small frequency segments in the computation. One key feature in positive 

20 species identification is the observance of a plurality of lines. If a line is missing, it 

generally implies that the species is not present. However, due to the strong interference 
in the spectrum by powerful neighboring lines, lines may either be previously marked as 
"present but distorted" or they may be altogether absent. Both cases are accounted for in 
the algorithm and the resulting spectrum is then fully diagnosed for the particular species, 

25 which if found also is assigned a concentration. 

The chemical species library used with the spectrometer is a pure database. It 
must be populated by spectral measurements nominally performed on the spectrometer, 
with care in preserving the calibrations. All the absorption lines are not necessarily 
stored in the library. In particular, weaker lines that are near to or in the noise (clutter) 

30 are excluded from the library because they would then confuse the recognition process 
that otherwise is amply accomplished by assessing only the species' predominant lines. 
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The library super-structure is adaptable depending on the applications. That is, the 
chemical order of importance is built into the hierarchy so that the most critical species 
for the application are examined first or more often. The library (preamble) and even the 
line content can be adjusted to meet these needs prior to use or in mid use, given suitable 
5 communications to the unit, or direct operator action if the unit is not in remote use. The 
subbands may also be used in the signal processing chain quite dynamically. Consider an 
assessment of a single subband that is quantifiably "poor." The system can then re-act by 
adding the next subband and refining its output. 

Fig. 1 1 shows the effect of clutter on performance of the spectrometer. Some 

10 applications of the spectrometer require performance in filed environments outside 
laboratory setting, in which false alarm avoidance is required. The use of the THz 
rotational domain is the largest single solution to this problem. That is, intrinsically, the 
clutter impediment is minimal for submillimeter spectroscopy to other device 
methodologies. However, clutter degrades the operation if a too limited line set of 

1 5 spectral data of chemical species stored in the library is used and/or if the net frequency 
interval used is inadequate. To a large extent, both can be controlled, and therefore do 
not represent a particular threat to the signal processing. The design concept includes a 
"tunability" feature that allows the system to be adjusted from higher detection 
probability to higher false alarm avoidance, as needed for the particular device 

20 application. It is noteworthy that the issue of false alarms should be put into the context 
of the application. The submillimeter spectrometer is intrinsically a very general purpose 
instrument. Given a very large library, a long term monitoring mission of a remote unit 
may simply issue a periodic report that is essentially a listing of atmospheric content. 
However, the spectrometer system may be set up to operate a "bell ringer" detector of 

25 some species or set of species. Only the first occurrence is relevant and false alarms must 
not be presented. So in this instance the system should be "tuned" away from the 
nominal knee of the ROC (nearest the upper left corner of the chart where the ideal 
system exists with PD=1 and PFA=0) downward toward the lower left region. All of this 
can be accomplished a priori or during the mission itself, affording the spectrometer 

30 systems the highest flexibility. 
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As described herein, the integrated spectral data can be combined with the GPS 
position/time and packaged for data communications. The values chosen are, for 
example, the averaged GPS coordinates and time during the corresponding measurement. 
Some application may require encryption for security purposes. The use of forward error 
5 correction (FEC) or other error correction codes (ECCs) may be used. It may be that a 
missed data point or two is non-damaging to the application and may not warrant the 
overhead and complexity of requiring error free communications. Some applications 
clearly can be best served by a collect and store methodology using a GB-sized disk 
module that is recovered after a sires of remote sensing. Encryption may still be 

1 0 performed in such applications. In the store and forward methodology, the packaged, 
integrated and registered (i.e., stamped with position, time and possibly other header 
information) data can be buffered for eventual transmittal to the A/C communications 
system for transmission to the remote processing site. In addition, physical parameters 
employed by the system (namely temperature and pressure readings, or volumetric values 

1 5 if variable volume techniques are required) that are required to calibrate the concentration 
of the gas sample must be returned with the sample data. The sample header format 
would extend to include such values for easy extraction and full time correlation with the 
sample. 

Other embodiments and uses of the invention will be apparent to those skilled in 
20 the art from consideration of the specification and practice of the invention disclosed 
herein. All references and written material cited herein for any reason are specifically 
and entirely incorporated by reference. It is intended that the specification and examples 
be considered exemplary only with the true scope and spirit of the invention indicated by 
the following claims. 
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